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Abstract 

The dependence of the subprocess cross section for dijet production at fixed trans- 
verse momentum on the (large) rapidity separation Ay of the dijets can be used to 
test for 'BFKL physics', i.e. the presence of higher-order (a s logAy) n contributions. 
Unfortunately in practice these subtle effects are masked by the additional, stronger 
dependence arising from the parton distributions. We propose a simple ratio test us- 
ing two different collider energies in which the parton distribution dependence largely 
cancels. Such a ratio does distinguish between the 'asymptotic' analytic BFKL and 
lowest-order QCD predictions. However, when subasymptotic effects from overall 
energy-momentum conservation are included, the BFKL predictions for present col- 
lider energies change qualitatively. In particular, the real gluon emission contributions 
to the BFKL cross section are suppressed. The proposed ratio therefore provides an 
interesting laboratory for studying the interplay of leading-logarithm and kinematic 
effects. 



Quantum Chromodynamics has been very successful in describing jet production in high 
energy collider experiments. Perturbative QCD at fixed order in the strong coupling constant 
a s provides sufficient predictive power for a wide variety of high energy phenomena. However 
in some regions of phase space, large logarithms can multiply the coupling, spoiling the good 
behavior of fixed-order expansions. In certain cases these large logarithms can be resummed, 
as in the Balitsky, Fadin, Kuraev and Lipatov (BFKL) equation ||, and the predictive 
power of the theory is then in principle restored. The latter statement should be subject to 
experimental test, since the predictions of BFKL are distinct from those of fixed-order QCD. 

Being asymptotic, the regions for which BFKL resummation is presumed to be necessary 
have until recently been beyond the reach of experiments, and therefore the predictions of 
BFKL have been difficult to test. The last few years, however, have seen experiments at the 
HERA ep collider and the Fermilab Tevatron that approach the kinematic regimes relevant 
to BFKL. At HERA the relevant regime is small parton momentum fraction x, where BFKL 
predicts a sharp rise in the structure function F 2 as x decreases. Further discussion can be 
found in Refs. 0, |3| and references therein; here we are concerned with BFKL physics at 
hadron colliders. 

At the Tevatron pp collider, and at hadron colliders in general, BFKL resummation 
applies to dijet production when the rapidity separation A of the two jets is large, A ^> 1, 
while fixed-order QCD should be adequate for A = 0(1). BFKL predicts, among other 
things, a rise in the dijet subprocess cross section a as A increases [§J], in contrast to the 
a — > constant behavior expected at lowest order. In practice such behavior can be difficult 
to observe because a gets folded in with parton distribution functions (pdfs), which decrease 
with A much more rapidly than the subprocess cross section increases. In particular, the 
BFKL rise with A gets killed by kinematic constraints. The challenge is to find measurable 
quantities in dijet production that are insensitive to the pdfs, but that retain the distinctive 
behavior characteristic of BFKL resummation. Some possibilities have been discussed in 
Refs. 0, ||, [7|, ||, H, [l(J, including the azimuthal decorrelation of the two jets: the multiple 
emission of soft gluons between the leading jets predicted by BFKL leads to a stronger 
decorrelation than does fixed-order QCD, and the prediction is relatively insensitive to the 
pdfs. 

Another possibility is to look for the increase in a with A by considering different collision 
energies Q. The idea is to choose A's that correspond to the same parton momentum frac- 
tions at different energies so that the pdf dependence is the same for both, thereby allowing 
the A dependence to be extracted. The D0 collaboration has in fact recently reported the 
results of a preliminary study of dijet cross section ratios at fixed parton momentum fraction 

in • 

In this paper we investigate the dependence of the dijet cross section on collision energies 
in leading-order QCD and in the BFKL approach. We define a cross section ratio at two 
collider energies in which the pdf dependence largely cancels, and in which the predictions 
of leading-order QCD and 'naive' BFKL can be distinguished. (By naive BFKL, we mean 
the predictions obtained by resumming leading-logarithm (a s A) n contributions arising from 
the emission of soft real and virtual gluon contributions in the absence of overall kinematic 
constraints.) However, one difficulty with the BFKL approach is that even for the largest 
A values which are accessible at present, subasymptotic effects, for example from energy- 
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momentum conservation, are likely to be important. This has led to 'improved' BFKL dijet 
cross section calculations, in which some of these effects are taken into account. In this 
study we will use the BFKL Monte Carlo approach developed in Ref. |10| (see also ||). We 
find that at present collider energies the predictions of naive BFKL are indeed substantially 
modified when subasymptotic kinematic effects are included. 

Our calculations will therefore give us a sequence of predictions — full leading-order QCD, 
asymptotic leading-order QCD, naive BFKL, improved BFKL — to confront experiment. 
In what follows we develop each of these approximations in turn. 

We begin with the lowest-order QCD inclusive two-jet production cross section as a 
function of the two jet rapidities yi,y2 arid their common transverse momentum p?- For 
simplicity we will consider the symmetric situation where the two rapidities are equal and 
opposite. Setting y 1 = —y 2 = \A gives for the differential cross section 



da 



dy 1 dy 2 dp 



i = 2fc^¥Ai ^ xf a (x,n 2 )xf b (x,n 2 )Y,\M(ab^cd)\' 
yi =- y2 =\& ^o07rp T cosn ^ 5 z\j aAc>d=q>g 



x = —j= cosn^Aj , 



(1) 

where J2 denotes the appropriate sums and averages over colors and spins and f a ,b{x, /x 2 ) are 
the parton densities. Here 

2pr 

and the subprocess matrix elements are functions only of the rapidity difference A. In the 
limit of large A the gg, qg and qq subprocess matrix elements squared become equal, up to 
overall color factors of C\, CaCf and C\ respectively. This defines the effective subprocess 
approximation to Eq. ([!]): 



da 



dy 1 dy 2 dp 2 T 



[xQ(x,fi 



2\12 



M2= i A 2567TJ4 cosh 4 (iA) 



E\M(gg^gg)\\ (3) 



where 

£ = S + ^E(9 + 9). (4) 

J q 

In an experiment (and in the BFKL formalism to be discussed below), we are interested 
in events with jets above some transverse momentum threshold Pt- Integrating over the jet 
transverse momentum pt > Pt then gives 



da 



dyidy 2 



E\M(gg^gg)\ 2 1 f^ 2 du 



m= _ m i A 256vr cosh 4 (|A) P| Ji u± 



xQ(x, /i 2 



x=Xu 

2' 



(5) 



with X = pT(min)/pr(niax) = 2Pt cosh(|A)/ v /i. In the asymptotic limit where A is large, 
the A dependence in the prefactor disappears so that 



E \M(gg 



' ~C 2 A c? s , (6) 



256tt cosh 4 (| A) 



;7T 
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and hence 



da 



<%id?/2 



yi=-y2- 



2P| 



x- 



du 2 



xQ(x, /i 2 ) 



x=Xu 



a F(X, fi 2 



(7) 



where 



2P 2 



and T contains the integration over the parton distribution functions. Equation (^) defines 
the asymptotic QCD LO cross section to which we will compare the BFKL predictions below. 

Fig. 1 shows the QCD LO dijet cross sections calculated using {%) the 2 — > 2 matrix 
elements (Eq. (H)), {%%) the effective subprocess approximation (Eq. (^)), and (Hi) the asymp- 
totic form of the effective subprocess approximation (Eq. (0)), for two pp collider energies 
y/s = 630, 1800 GeV. The jet p T threshold is P T = 20 GeV/c, the partons are CTEQ4L 
with p, = P T |I2]| , and a s is evaluated in leading order also at scale fi = P? (the numerical 
value is 0.170). The effective subprocess approximation is seen to reproduce the exact result 
to within a few per cent over the entire A range. The asymptotic form is approached at 
large A, as expected. 

At higher orders in QCD perturbation theory the subprocess cross section receives large 
logarithmic corrections from soft gluon emission in the rapidity interval between the two 
leading jets. For fixed coupling a s we have 



a = cr 



1 + an (a s A) n + 



n>l 



(9) 



with o" as defined above. The BFKL formalism to be discussed below resums the logarithms 
and in the leading-logarithm, fixed-coupling 'naive' approximation gives a subprocess cross 
section which is not constant but grows asymptotically with A: b ~ exp(AA) with A = 
4C A ln2a s /7r w 0.5. 

At fixed y/s and minimum transverse momentum Pt, both a and T depend on A. In 
fact because of the shape of the parton distributions the latter quantity decreases rapidly 
with increasing A and vanishes at the kinematic limit A = 2 coshT 1 (\/s / 2Pt) ■, as in Fig. 1. 
It is therefore difficult to observe the relatively slow rise with A of a (see for example Fig. 8 
of Ref . ) . The original idea of Ref . [Q was to increase the collider energy y/s as A increases 
such that X and therefore JF remains fixed. Any observed rise in the cross section could 
then only arise from higher-order contributions to a. 

If two collider energies are available [\/s 1 and y/s 2 say) one can make use of this idea 

by comparing the dijet cross sections at two rapidity separations (Ai and A 2 ) for which the 

asymptotic leading- order cross sections are equal. Specifically, if for a given Ai we define A 2 

such that _ 

cosh(|Ai) = ^ 

cosh(iA 2 ) y/s 2 

then the cross sections defined by Eq. (0) will be equal, i.e. da(y/s 1 ,A 1 ) = da(y/s 2 , A 2 ). 
Using for y/s~i and y/s~2, respectively, the two Tevatron collider energies 630 and 1800 GeV, 
we have the following values for Ai and A 2 : 
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A 1 ( v /I7 = 630 GeV) 


0.0 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


6.8 


A 2 ( v /^= 1800 GeV) 


3.42 


3.68 


4.33 


5.19 


6.13 


7.11 


8.10 


8.90 



Note that for large Ai Eq. fllCf ) reduces to A 2 = Ai + ln(s 2 /,si) = Ai + 2.1, as can be seen 
in the table. 

The equality of the asymptotic cross sections at lowest order for the two A's leads us to 
define a cross section ratio 

da(^ v Ai) 



R 



12 



A, 
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do-(Vs 2 »' 

as a function of Ai, with A 2 given by Eq. (|1(J). This ratio has the advantage of being directly 
accessible experimentally, since it depends on the rapidity difference between the two 'out- 
side' jets, a quantity more easily measured than for example the parton momentum fractions 
themselves. By construction, i? 12 = 1 in asymptotic LO QCD, but notice that the ratio is 
not equal to unity for the leading-order cross section at subasymptotic rapidity separations 
calculated using exact matrix elements, or using the effective subprocess approximation. 
This is because the prefactor multiplying the integral in for example Eq. (|5|) only becomes 
independent of A at large A. This is illustrated in Fig. 2, which shows i?i 2 as a function of 



Ai for y /s 1 = 630 GeV and y/s 2 = 1800 GeV, using the parameters of Fig. 1. The effective 
subprocess approximation curve is very close to the exact curve and is not shown. The ratio 
approaches unity at large A 1; as expected. Not surprisingly, the 'exact matrix element' ratio 
is insensitive to the choice of parton distributions. The dashed line shows the ratio obtained 
using GRV94LO distributions 



m. 



The curves are almost identical. Similarly, the scale 
choice dependence is also very weak. For example, using the more natural 'running' choice 
fi = Pt instead of fi = in the parton distributions and in a s has a negligible effect on the 
ratios. 

The BFKL formalism resums the leading-logarithmQ (a s A) n higher-order corrections 
which arise from multiple emission of real and virtual gluons in the rapidity interval be- 
tween the two leading jets. The theoretical details of how this is achieved can be found 
elsewhere (see for example Ref. ||15|| ) and will not be repeated here. 

In fact the naive BFKL prediction, in which for example subasympotic effects from overall 
energy momentum conservation are ignored, has an analytic representation which leads to a 
simple expression for i? 12 in the large Ai limit. If the coupling a s is assumed to be constant, 
all kinematic constraints are ignored, and only the leading logarithms are resummed, then 
the result is 

fa s C A 

O"0 L'O 



a 



71 



12 



with 



C (t) 



Xo{z) 



1 

2^ 
Re[V(l) 



dz 



z 2 + \ 



IZ) 



(13) 



1 In fact, sub-leading logarithms (a™A n 1 etc.) can also be resummed in principle (l4|. In practice, 
however, only the leading contributions are known at present in a form which is phenomenologically useful. 
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The asymptotic (large t, equivalently large A) behavior is 



c ° (t) ~ T^r 4tog2 ' ■ (14) 



It follows immediately that the naive BFKL prediction for the ratio Ru at large rapidity 
separations is 

Rl2 = °) ; . (is) 

C (^A 2 ) 

The prediction is shown in Fig. 2. Note that now R12 is well below the QCD LO curves, 
and in fact at first sight this appears to be a primary test of the presence of higher-order 
'BFKL-like' corrections to the dijet cross section. The formalism is of course only supposed 
to be valid for large Ai (which implies large A2 also) and so we only show the predictions 
for Ai > 2, which is roughly where the leading-order prediction begins to approach its 
asymptotic limit R12 = 1. The asymptotic large Ai behavior of the ratio is readily obtained 
from Eqs. and (|15|): 

Ru^(^)\ A=^41og2. (16) 

One interesting aspect of this result is that one is tempted to use the experimentally mea- 
surable quantity /C = log(i?i2)/log(si/s 2 ) to determine the effective BFKL exponent A, 
i.e. 

^-IlO asymp. = ^ ' 
KlnaiveBFKL - ^41 Og 2 = 0.45, (17) 

where the numerical value is obtained from using a s evaluated at /x = Pt = 20 GeV, as in 
Fig. 2. However it is clear from the figure that the asymptotic behavior is attained only very 
slowly. For this choice of parameters we expect R\2 — > 0.39 as Ai — ► oo, but we see that the 
ratio has only decreased to R\2 = 0.45 at the kinematic limit Ai = 6.9. Furthermore, as we 
shall see below, an improved BFKL calculation does not lead to such a simple prediction. 

So far we have considered the predictions of naive BFKL, which can be written ana- 
lytically but at the cost of making several assumptions which do not hold in experimental 
situations. In particular, the analytic solutions result from integrating over arbitrarily large 
values of the transverse momenta of emitted gluons, and the analytic phase space integra- 
tions at the subprocess level prevents a proper incorporation of parton distributions. In 
addition, a s is assumed to be fixed. The BFKL approach can be improved and these as- 
sumptions avoided by recasting the BFKL cross section as an event generator [|, [TOj] using 



an iterated solution to the BFKL equation. Kinematic constraints and the running of a s 
can then be included. The effect is to restrict the growth of a with increasing A compared 
to naive BFKL, since the emission of relatively large transverse momentum gluons, which 
give a positive contribution to the resummed cross section, is suppressed. Details of the 
calculation and application to the azimuthal decorrelation can be found in [|l(| . 
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Here we repeat the calculation of the cross section and ratios Ru using the Monte Carlo 
of Ref. [0] to get an 'improved' BFKL prediction. One feature of the MC calculation 



that will be relevant to our results is worth mentioning. The MC solution to the BFKL 
equation is obtained by separating the contributions from real gluon emission into 'resolved' 
and 'unresolved' contributions, the idea being that below a certain energy scale /i , emitted 
real gluons are not detectable in practice. The contribution from unresolved real gluons is 
combined with the virtual gluon contribution to give an overall suppressing form factor. The 
differential cross section takes the form, for fixed 



da gg _ a 2 C\ 



Pti 



d 2 p T id 2 pT2dAy Vt\Vt2 
and TZ can be written as a sum over resolved real gluon emissions: 



o s C,4A/7r 

n{p T1 ,pr 2 ,Ay) , (18) 



TZ(p T1 ,p T2 ,A) = ]T7^(p T1 ,p T2 ,A) , (19) 



n=0 



with 7Z^ = 5(pti + PT2) and TZ^ for n > given in [TO]. The point is that even for n = 
the form factor t 2 - appears in the cross section and represents a suppression of 

Pti] 

the probability that no resolvable gluons will be emitted. The inclusion of resolvable gluon 
contributions (n > 1) counteracts this suppression and removes the dependence on fiQ. 

Fig. |H shows the dijet cross section as a function of A for the naive BFKL and improved 
BFKL MC cases at the two collision energies, with Pt = 20 GeV. Asymptotic QCD LO is 
also shown for reference. There are several features worth noting here. First, the naive BFKL 
cross section (dashed curve) is always largest, because it includes the analytic subprocess 
cross section given in Eq. (|l|) , which allows emission of any number of gluons with arbitrarily 
large energies. The curve falls off rather than increasing because a is multiplied by parton 
densities, but even those incorporate only lowest order kinematics in this case. 

When exact kinematics for entire events are included in both the subprocess cross section 
and the parton densities, as in the BFKL MC (solid curve), there is a dramatic suppression 
of the total cross section.0 In fact the suppression is so strong that it drives the BFKL MC 
cross section below thai for asymptotic QCD LO. The reason is due to simple kinematics: the 
QCD LO cross section contains only two jets, but the BFKL MC cross section also includes 
additional jets, each of which increases the subprocess center-of-mass energy and elicits a 
corresponding price in parton densities. In the naive BFKL calculation, the contribution 
to the subprocess energy from additional jets is ignored and their net effect is to combine 
with the virtual gluons to increase the subprocess cross section. Interestingly, the effect 
is compounded by the fact that the (naive) BFKL increase of a only starts at 0(a 2 ) in 
perturbation theory, i.e. Co = 1 + 0(a 2 ) in Eq. (0). The 0(a]) correction is zero because 
the real and virtual gluon contributions exactly cancel. However the effect of the pdfs is to 
suppress only the real contributions, giving an overall negative correction at this order. 



2 We show the fixed a s expression for simplicity; see JT(J for the corresponding running a s expression. 
3 The running of a s , which we include, also contributes to the suppression, but it has a much smaller 
effect than kinematics. 
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The ratio R\ 2 as calculated in the improved BFKL MC is shown in Fig. |], again with the 
naive BFKL and asymptotic QCD LO predictions. Here the effects of kinematic suppression 
and the parton distributions are quite dramatic. The BFKL MC curve does not fall between 
the naive BFKL and QCD LO curves as one might expect. Instead, it is everywhere greater 
than or equal to the QCD LO curve, and in particular, i? 12 is everywhere greater than 1. 

This behavior can be understood as follows. We can represent Ri 2 schematically as 

Ru ~ 1 + + , (20) 

1 + a s A 2 C 2 + . . . ' v ; 

where the Ci are coefficients modified by the effects of the kinematic suppression. By the 
arguments given above, we expect Ci < 0. In the formal^ limit Ai — > 0, the nonleading 
terms in the numerator vanish, the denominator is < 1 and hence Ri 2 > 1. On the other 
hand, in the limit Ai — > A max , there is no phase space available for any gluon emission, 
hence Ci — > and R12 — > 1. The observation of a 'BFKL effect' would require a region of 
phase space where A» ^> 1 but Ci ~ 1. Increasing the collider energy at fixed transverse 
momentum threshold would eventually lead to this situation. In such a region we would 
expect to see R\ 2 < 1, as predicted by the naive BFKL calculation shown in Fig. £|. Finally, 
why does the BFKL MC curve in Fig. |] rise again at large A 1; instead of approaching 1 as 
simple kinematics would suggest? The origin of this behavior arises in the form factor (see 
Eq. (|T8D) which suppresses the probability of no real (i.e. resolved) gluon emission in the 
rapidity interval between the dijets. As A increases towards its maximum allowed value there 
is an effective upper (kinematic) limit on the transverse momentum of each emitted gluon. 
This in turn leads to a suppressing form factor with the 'artificial' parameter /x replaced by 
a physical parameter determined by kinematics and the pdfs. In other words, the all-orders 
leading-log contribution from virtual gluon emission that is built into the overall form factor 
is not completely compensated by corresponding real emission because the kinematics do 
not allow it. The kinematic suppression increases with A, hence A 2 > Ai implies R\ 2 > 1. 

The arguments above illustrate the general behavior of the cross section ratio - starting 
out at a value greater than 1, falling and then increasing again at large A - as predicted by 
the 'improved' BFKL Monte Carlo calculation. The details, however - how far the ratio falls 
(and in particular whether it falls below unity as predicted in the naive calculation) and for 
how long - depend on the specifics of the experimental configuration, such as the collider 
energies and the transverse momentum threshold. In this study we have restricted ourselves 
to parameter values which are currently accessible at the Tevatron. Clearly it would be 
desirable to have more than the two Tevatron pp energies available for such a measurement. 
In principle, the LHC at 14 TeV and RHIC in pp mode at 500 GeV would expand the range 
of A. In practice, one is limited by A max at the lower energy. Furthermore, comparing jet 
measurements at different machines with different detectors can be challenging, and ideally 
one would like to construct ratios of measurements made at a single machine in a single 
detector, so that some systematic uncertainties would cancel. For example, running the 
LHC at 10 and 14 TeV would make such a measurement possible, although ideally the 
two energies should be more different. The essential point is that higher collider energies 

4 We do not of course expect the BFKL formalism to apply for small rapidity separations; here our aim 
is simply to explain the behavior of the predictions of the BFKL MC calculation in Fig. ||. 
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appear to be necessary in order to prevent the BFKL effects being swamped by kinematical 
constraints. 

In summary, we have shown that the effects of the increase in the dijet subprocess cross 
section predicted in the naive BFKL approach can in principle be detected by measuring 
the ratio of cross sections at energies and rapidities chosen such that the asymptotic QCD 
LO ratio is equal to 1. Forming the ratio minimizes the effects of the parton distribution 
functions. However an improved BFKL calculation which includes the subasymptotic effects 
of the conservation of energy and momentum gives a result for the ratio that is qualitatively 
different from that of naive BFKL. In particular, for dijet production at the two Tevatron 
energies, the improved BFKL calculation gives either agreement with lowest-order QCD or 
deviations opposite to those predicted by naive BFKL. If our arguments about the kinematic 
suppression of higher-order real gluon emission are correct, then it would appear that a fixed- 
order perturbation theory approach would be a more appropriate calculational framework for 
the cross section ratio at present collider energies than all-orders resummation. Therefore, 
an important next step is to perform a calculation of R12 in exact next-to-leading-order 
QCD, to determine whether it provides a satisfactory description of the data. 

Finally we emphasise that this study has concentrated on the cross section ratio only. 
It may well be that other quantities are better suited to revealing underlying BFKL effects. 
Indeed, we note that the azimuthal decorrelation of the dijet pair || |], as calculated in 
the naive BFKL approach, is already non-zero at 0(a s ) 0, and that its effects are not 
completely swamped by kinematics in an improved BFKL calculation [pL0| . 
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Figure 1: The dependence of the leading-order (2 — > 2) dijet cross sections on the dijet 
rapidity separation. The partons are the CTEQ4L set ||12|| . The three curves at each collider 
energy use: (i) exact matrix elements (solid lines), (ii) the effective subprocess approximation 
(dashed lines), and (iii) the asymptotic (A ^> 1) form of the latter (dotted lines). 
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1.2 




; naive BFKL p T >20GeV/c [ 

i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

1 2 3 4 5 6 7 



Figure 2: The ratio Ri 2 of the dijet cross sections at the two collider energies = 630 GeV 
and y/s 2 = 1800 GeV, as defined in the text. The curves are: (i) the exact leading-order 
(2 — > 2) predictions using CTEQ4L (solid curve) |12| and GRV94LO partons (dashed curve) 
jl3|1 , both with /i = P T = 20 GeV, and (ii) the 'naive BFKL' prediction (dot-dashed curve). 
Note that the asymptotic leading-order prediction is Ri 2 = 1. 
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Figure 3: The dependence of the BFKL and asymptotic QCD leading-order dijet cross 
sections on the dijet rapidity separation. The pdfs are the CTEQ4L set ||12|| . The three 
curves at each collider energy use: (i) 'improved' BFKL MC (solid lines), (ii) 'naive' BFKL 
(dashed lines), and (iii) the asymptotic (A ^> 1) form of QCD leading order (dotted lines). 
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Figure 4: The ratio R\2 of the dijet cross sections at the two collider energies y / i 1 = 630 GeV 
and y / i 2 = 1800 GeV, as defined in the text. The curves are: (i) the 'improved' BFKL MC 
predictions using CTEQ4L |T2j pdfs (solid curve), with fi = Pt = 20 GeV, (ii) the 'naive' 
BFKL prediction (dashed curve), and (iii) the asymptotic leading-order prediction (dotted 
curve) Ru = 1. 
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